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^s^ ■ Abstract: This paper continues the study of the Gregory-Laflamme instabihty of black 

strings, or more precisely of the order of the transition, being either first or second or- 
der, and the critical dimension which separates the two cases. First, we describe a novel 
Oh. method based on the Landau-Ginzburg perspective for the thermodynamics that some- 

r~| , what improves the existing techniques. Second, we generalize the computation from a 

circle compactification to an arbitrary torus compactifications. We explain that the criti- 
cal dimension cannot be lowered in this way, and moreover in all cases studied the transition 



\^ ' order depends only on the number of extended dimensions. We discuss the richer phase 

structure that appears in the torus case. 
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1. Introduction and summary 

In the presence of a compact dimension Gregory and Laflamme (GL) discovered that uni- 
form black strings are perturbatively unstable below a certain critical dimensionless mass 
density ||l[ . The order of the transition can be computed by following perturbatively the 
branch of non-uniform solutions which emanates from the critical GL string, as first shown 
by Gubser in the case of a five-dimensional spacetime Q where the transition is first order. 
That calculation was generalized by one of us (ES) to arbitrary spacetime dimensions with 
the surprising result that the transition is first order only for D < D* = "13.5" while 
it is second order for higher dimensions Q. Here first order means a transition between 
two distinct configurations, while a second order transition is smooth - the uniform string 
changes smoothly into a slightly non-uniform string. Kudoh and Miyamoto repeated 
the calculation in the economical Harmark-Obers coordinates [^], confirmed previous re- 
sults and observed that in the canonical ensemble the critical dimension actually changes 
from D* = "13.5" to D*^^ = "12.5". All this data constitutes an important piece in the 
construction of the phase diagram for this system (see |^ and Q for a review) . 



The present paper includes two main results. First, we show how to somewhat improve 
the existing method of calculating the transition order by employing a Landau-Ginzburg 
perspective (the basic idea was described already in Appendix A of [^). Secondly, we 
generalize from the usual S^ = T^ compactification to an arbitrary torus compactification 

Landau-Ginzburg improvement to the method. In Section g we review our set-up and 
how the description of phase transitions is achieved by the Landau-Ginzburg (LG) theory, 
where one expands the free energy of the system around the critical point in powers of 
order parameters. In particular, it is known that as long as the coefficient of a certain 
cubic term in the free energy is non- vanishing then the transition is first order. If the cubic 
term vanishes, for instance due to a parity symmetry such as in our case, then it is the sign 
of the coefficient of a certain quartic term, which we denote by C, that determines whether 
the transition is first order or higher (of course if this term vanishes one has to go to higher 
terms). A simpler and more intuitive form of this criterion, which is, alas, less general and 
does not apply in our case (as expanded in the text), is that the transition is second order 
iff the critical solution is a (local) minimum of the free energy. 

Before we can compare the Landau-Ginzburg method with Gubser's method, we should 
recall the features of the latter. There one computes order by order the metric of the static 
non-uniform string branch emanating from the critical GL string. The first order is nothing 
but the GL mode. At the second order one computes the back-reaction. Finally the third 
order is computed, or more precisely only the first harmonic along the compact dimension, 
from which one can finally compute the leading coefficients of the changes in mass and 
entropy, r]i, a"2 of the new branch. The sign of these two quantities is correlated and 
determines the order of the transition. 

At first sight the two methods look quite different. However, we show that in the LG 
method one also needs to precisely compute the second order back-reaction to the metric. 
The third order however is not required in LG (thereby avoiding the solution of a set of 
linear differential equations with sources). Rather one needs only to expand the action to 
an appropriate quartic order, to substitute in the results from the first and second orders 
and perform certain integrals that add up to the constant whose sign determines the order. 

A way to understand the simplification is the following: in Gubser's method one com- 
putes the third order, but it turns out that all that is really needed is the projection of 
the third order onto the GL mode. That is precisely the reason why the first harmonic 
sufficed (as the GL mode is in the first harmonic). Our substitution into the quartic order 
of the free-energy achieves exactly that, without the need to compute other properties of 
the third order. 

In Section ^ we perform the "Landau-Ginzburg" calculation for an S^ compactification 
in various dimensions and verify that we get the same bottom- line coefficients ryi, cr2 as in 
the previous method, see table ^ 

Torus compactification. It is interesting to generalize the compactifying manifold, and the 
simplest option beyond the circle S^ = T^, is a product of circles, or more generally a 



p-torus T^. The number of extended spacetime dimensions is denoted by d and the total 
spacetime dimension is D = d+p. The critical GL density for such a torus compactification 
is easily found to be given in terms of the shortest vector in the reciprocal lattice |^. 

We proceed to analyze the transition order in Section^. First, we motivate restricting 
ourselves to square torii. Basically, we view the space of torii as having two boundaries - 
on the one hand highly asymmetrical torii, where one (or more) dimensions are much larger 
than the rest, and on the other hand highly symmetrical torii such as the square torus. 
Since the limit of a highly asymmetrical torus reduces to the case of a lower dimensional 
torus (mostly the well- understood case of §^ compactification) , we argue that by studying 
the opposite limit of a highly symmetrical torus, we achieve an understanding of both limits 
and thereby also some understanding of the intermediate region of general torii. 

For a square T^ torus compactification, p modes turn marginally tachyonic at the 
same (GL) point. We find that the constant C is replaced hy a p x p quadratic form C*-', 
in order to allow for the various possible directions in the (marginally) "tachyon space", 
and that the transition is second order iff C"^ is positive for all directions. Namely, it is 
enough that there is a single direction in tachyon space which sees a first order transition 
for the transition to be one. Taking into account the T^ results we may immediately deduce 
that the critical dimension cannot be lower than in the T^ case with the same number of 
extended dimensions, d. This shows that the indicators of [^ for second order transition 
at lower dimensions, which were part of the motivation for this work, were misleading, as 
further discussed in the text. 

Due to the high degree of symmetry of the square torii C*-^ consists only of two indepen- 
dent entries: all the diagonal entries are the same, denoted by C^ , and all the off-diagonal 
entries are the same, denoted by C^. Since the diagonal term is precisely the one computed 
in the T^ compactification, we set to compute the off-diagonal term. Due to the symmetry 
C^ is the same for all p and for that purpose it suffices to consider p = 2, namely we 
consider the square T^ torus. The only parameter remaining is the number of extended 
dimensions. 

In practice we do not compute C^ directly, but rather the effective C for a diagonal 
direction in tachyon space, which we denote by C. Then we solve for C^ which is a linear 
combination of C^, C\p=2- 

In Section || we proceed to the actual calculation for T^. Once we have chosen the 
diagonal direction in tachyon space we are not bothered any longer by the presence of several 
(marginally) tachyonic modes. However, the number of metric components involved in the 
calculation (back reaction and quartic coefficient) is larger than in the T^ case. Certain 
discrete symmetries are found to be helpful in simplifying the calculation. The calculations 
and results are described in detail, thereby applying our improved method to achieve new 
results. 

In Section |6| we analyze the results and their implications. We find that for all the 
studied values of d where the T^ transition is second order, the T^ transition is also second 
order for all p. Combining this result with observation mentioned above we conclude that 
the transition order for square torii shows some robustness in that it depends only on d, 
the number of extended dimensions, and not on p. 



In addition we discuss some subtler implications, including the finding that for almost 
all d the diagonal direction in tachyon space is disfavored relative to turning on a tachyon 
in a single compact dimension, and in this sense we have spontaneous symmetry breaking. 

The appendices |A|, ^ contain the details of the full equations, of the numerical tech- 
nique and summarize various numerical parameters. 

2. Set-up and Landau-Ginzburg theory of phase transitions 

2.1 Set-up 

The Gregory-Laflamme instability and associated phase transition physics appear in the 
presence of compact dimensions, namely for backgrounds of the form M '^ x Y^ where 
yP is any p-dimensional compact Ricci-flat manifold, d > 4 is the number of extended 
spacetime dimensions, and the total spacetime dimension is D = d + p. 

The simplest compactifying manifold is a single compact dimension y = S^. We shall 
first demonstrate the Landau-Ginzburg method in that context and then turn to y = T^, 
a square p-dimensional torus. We denote the torus size by L and work with Euclidean time 
of period (3 (corresponding to a canonical ensemble). Such backgrounds are characterized 

by a single dimensionless constant 

(3 , , 

fi^oc- , (2.1) 

where henceforth we shall omit the subscript (3. 

The non-rotating black objects in which we are interested are static and spherically 
symmetric (in the extended dimensions). After suppressing the time and the angular co- 
ordinates the remaining essential coordinates are r, the radial coordinate in the extended 
spatial directions, and z*, 1 < i < p the periodic coordinates z* ~ z* + L which parame- 
terize T^. Thus, the essential geometry has p + 1 dimensions and we employ "cylindrical 
coordinates" (r, 2;*). 

The uniform black p-brane^, the background around which we perturb (see figure |l]) 
is given by 

d'' = ^^Schw + ^'Y (2.2) 

where dsy is the metric on Y^, which for our square torus is 

p 
dsl = ^dz'dz\ (2.3) 

and ds% , is the d-dimensional Schwarzschild black hole, 
bchw ' 

1 



^^Schw = -/(^) ^* +T7^^d^ + ^ d^d-2, (2-4) 



where 



/(r) = 1 - f'^V"' , (2.5) 



^A 1-brane is often called "a string" 



E^ 



Figure 1: The uniform black string together with the definition of the cylindrical coordinates (r, z). 
To is the Schwarzschild radius. 
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and c?r2^_2 is the metric on the unit sphere, rg is related to the black hole mass, M, via 

m 

(2.6) 
1:7-7- is the area 






" (d - 2) nd-2 ' 

where Gd is the d-dimensional Newton constant, and ^d~i — "Trf/W 

of a unit sphere S . In Euclidean signature the space-time ends at the horizon. As 
usual, requiring the absence of a conical singularity there fixes (3, the asymptotic size of 
the Euclidean time direction, and it is given by 



r(f) 



27r 
P = — 
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(2.7) 



Gregory and Laflamme discovered that the uniform black string solution ( |2.4D de- 
velops a z-dependent metric-instability below a certain critical (dimensionless) mass |[l|. 
Equivalently, for a fixed tq there is a critical wavenumber 



k = 2tt/L 



(2.8) 



for instability. At A; = kcL the GL mode is marginally tachyonic, namely a zero-mode. 
In p], H] the critical GL lengths were obtained for Schwarzschild black holes in various 
dimensions (see table ||). From these the high d asymptotic form was extracted and later 
proven analytically in Q to be 

kcL^Vd—. (2.9) 

This means that for large d the black string becomes unstable when r^/L ~ ydj namely 
when it is quite "fat" and this indicates that such a string would not decay into a black 
hole which would not "fit" inside the extra dimension. 



d 


4 


5 


6 


7 


8 


9 


10 


11 


kcL 


.876 


1.27 


1.58 


1.85 


2.09 


2.30 


2.50 


2.69 


d 


12 


13 


14 


15 


19 


29 


49 


99 


kcL 


2.87 


3.03 


3.19 


3.34 


3.89 


5.06 


6.72 


9.75 



Table 1: The marginally static mode wavenumbers kcL in units of r^ as a function of d [^ ^ . 



Having an unstable mode brings us to the issue of the order of the transition and the 
existence of additional phases. Rather than review the state of the art in this respect, we 
proceed in the next subsection to describe a somewhat novel method, which in section ^ 
will be shown to reproduce the known results and then will be used in section |5| to obtain 
new results. 

2.2 Landau-Ginzburg: review and application 

The central insight of the Landau-Ginzburg (LG) theory of phase transitions (see for exam- 
ple [|l^]) is that the nature of a phase transition can be deduced from the local behavior of 
the free energy around the critical solution. This local analysis is achieved by focusing on 
the low energy modes and zooming, namely carrying a Taylor expansion up to a specified 
order. For pedagogical reasons^ we divide the discussion into three steps: first we discuss a 
one-dimensional configuration space which is the simplest example, then the general con- 
figuration space and finally we utilize the translation symmetry in z of our background to 
arrive at the final form of our formulae. 

One- dimensional configuration space. The main idea can be demonstrated by a system 
with a configuration space consisting of a single variable A and by a control parameter /x. 
The thermodynamics is encoded by the free energy function F = F(X; /i). Generically the 
free energy does not have a definite parity in A, however sometimes additional symmetries 
of the system make the free energy even in A, namely F{—X; //) = -F(A; //), which will be 
seen to hold in our case. This reflection symmetry implies that A = is an extremum of 
the free energy (namely that the linear term in the Taylor expansion of -F in A around 
A = vanishes). Assuming the existence of a critical solution with a marginally stable 
mode means that for some critical value of // denoted by /ic the quadratic term vanishes. 
Thus F{X; ^) = Fo{iJ,) + A{^ — fic) A^ -|- . . . where A is some constant, which we take to 
be positive without loss of generality (if A were negative we could redefine fj, — > — /i). In 
order to determine the order of the transition it is sufficient to expand F further as follows 

F(A; fi) = Foifi) + A6nX^ + C X'^ + . . . (2.10) 

where we defined 

5fi = fj.-Hc (2-11) 

and C is another constant. Actually the whole function -Fo(m)) is not required for the 
LG theory, only its expansion around /u ~ ^c) but we shall not be concerned with this 
expansion. 

Once the constants A, C are determined the local thermodynamics may be deduced 
by following the form of the free energy and its extrema as /u is adjusted in the vicinity of 



/ic- Seeking additional extrema of (2.10) through dF/dX = results in a new branch of 
extrema in addition to A = namely, 



-2 A{^i-Hc 



xi = - \c ' ^^-^^^ 



^A brief outline of the method was already present in appendix A of 
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Figure 2: An illustration of a first order phase transition. A condensed plot showing two kinds of 
graphs. The thin solid lines show the free energy as a function of A for a sequence of fi values, while 
for the thick lines the vertical axis becomes /i (the horizontal remains A) and they designate the 
various phases corresponding to the extrenium of the free energy. The free energy has a minimum 
for fi > ^c that corresponds to the stable symmetric phase (thick solid line) which becomes unstable 
below fic (thick dashed line) . It follows that the asymmetric phase branch emergent from the critical 
point (thick dotted line) is unstable since the free energy has a negative direction for /i > /ic at 
A = Afc. Note that the free energy in this phase is higher relative to that of the critical state. 



where A^ = Af,(^) describes the new branch. Note that since A^ > the new branch exists 
only for sgn(/i — fic) = — sgn(C) (locally). Depending on the sign of C we now get two 
possibilities which are depicted in figures ^, |3[ For negative C the new branch exists only 
for fj. > fic and thus cannot serve as the end-state of the decay for fi < fic (moreover it is 
unstable) , and therefore a first order transition will occur into a phase at finite distance in 
configuration space. For positive C the new branch can and does serve as the end-point of 
the decay. The transition is smooth as one lowers fj, below fic- The difference in free-energies 
between the two phases is given by 



Fbit^) - Foifi) = F(Afe(/i);^) - Fo(m) 



A^ (/i - fic 
AC 



-C\^ 



b ■ 



(2.13) 



In particular there is a jump in the second derivative of F^ — Fq with respect to fi at fi = fic- 
Hence by definition this transition is second order. 

Given A, C and hence the form oi F = F{fj) , which is the fundamental thermodynamic 
potential (in the canonical ensemble) we may find the rest of the thermodynamic quantities 
through the usual thermodynamic relations, as we now describe. 

We particularize to the black branes in which case the temperature is T = h/ (3, and so 
the dimensionless temperature is ^ = fi~^ ( p.l[ ). Rephrasing ( p. 12 ) we obtain the leading 
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Figure 3: An illustration of a second order transition from a symmetric to an asymmetric phase, 
with the same condensed conventions as in figure 0. The free energy (designated by the thin solid 
line) has a minimum for ii > iic corresponding to the stable symmetric phase (thick solid line) which 
becomes unstable below iic (thick dashed line). The free energy is at minimum also for fi ^ jic and 
the minimum continuously moves away from A = 0, indicating that the emergent asymmetric phase 
is stable and has the asymmetry developing smoothly. 



behavior for the dimensionless temperature of the emergent branch 

6e__ 5jj^_ 20 
6 fl A fiGL ' 






hxl 



(2.14) 



The entropy and the mass of the solutions are S = fi'^d^F and M = d^{p.F). We can 
find the coefficients of the second order variations of these quantities relative to the critical 
phase 






si = 0{X^) of fi^d 



mi = 0(a2) of - d, 



f,[F{^i)-Fo{fi)]\ _iiglA 



(2.15) 
(2.16) 



^ F{ijl) J d-2 

The variation of the dimensionless mass, 

ri = GdM/L'^-^ = GdMfi'^-^/P'^-^ 

in a canonical ensemble for given (3 is comprised of two terms arising from variation of M 
and n 



(2.17) 



rii = ( 0(A2) of ^ J = -{d - 3)9i + mi 



(2.18) 



Finally, we can also determine the entropy difference between non-uniform and the 
uniform phases with the same mass P, ^ 

'-'non— uniform -. , \ 2 , \ 4 i 

— = 1 + aiXf, + a2Xb + ■■■, 

"-'uniform 

d-3 d-2 

Note that vanishing of ai is guaranteed by the first law. 

Higher dimensional configuration space. While the discussion above reveals the essential 
point, it considers the free energy to be a function of a single configuration variable, while 
in our case it should be considered a function over the space of all metrics with fixed 
boundary conditions, which is an infinite dimensional space. More explicitly 

-PF[g^,]=lGH[g,u]^j^y^R + 2j^jK-K'']^ , (2.20) 

where the first integral is the bulk contribution of the Ricci scalar and the second is a 
boundary contribution where K is the trace of the second fundamental form on the bound- 
ary, and Kq is the same quantity for a reference geometry. The space of metrics considered 
is that of all metrics which asymptote to the reference geometry, which in our case is flat 
R X Sl X T^ characterized by dimensions d and p and the asymptotic constants /3, L. 
We denote by z^ the coordinates along the compact dimensions, and by r the radial 
coordinate. The most general metric which is static and spherically symmetric (in the 
extended dimensions) 

ds^ = e^^dt^ + ds^^^^r) + e^^dnl_^ , (2.21) 

where A and C are functions of (r, z*), ds"^ i^ is an arbitrary metric on the (r, z*) space 
and since the metrics are static we might as well work with Euclidean signature. 

We denote the metric fields collectively by x. We are interested in a background which 
is the critical GL string, denoted by x^^\ and we decompose the metric fields into the 
background, x^^\ and fluctuations, denoted by^ X 

x = x^^^+X. (2.22) 

Since we wish to follow the new branch emanating from the zeroth order solution, x^^\ 
we expand in a perturbative parameter A 

oo 

X = Y^ VX» = A X(i) + A^ X(2) + . . . (2.23) 

The first order contribution must be a zero mode, namely the GL mode Xql 

X(i) = Xgl , (2.24) 

^In general X can be considered to be a vector. While in our case it is infinite dimensional (since the 
fiuctuations, X, are a set of metric functions, and as such represent infinitely many modes), more generally, 
if the configuration space is finite-dimensional then X would be a finite dimensional vector X*. 



and it will be sometimes convenient to denote the second order term by Xbr, the back- 
reaction 

X(2) = Xbr . (2.25) 

We need to expand the free energy F to fourth order in A. To that purpose we first 
expand F in X 



FiX)=Fo + F2iX,X)+F3iX,X,X)+F^iX,X,X,X)+... 
, OF 



(2.26) 



where Fi are multi-linear expressions^. The linear term vanishes since by assumption X = 
is a solution. 

Next, we decompose the fluctuations^ X into the marginal GL mode A, and all the 
rest, X 



X = (A,X) 

X =~XXgl + X 



(2.27) 



In terms of the perturbations series ( 2.23| ) we have 

A = A 
X = X Xbr + ■ 



(2.28) 



So A and A are equal along the perturbation path, and the separate notation is intended 
to distinguish between the direction in the space of metric fluctuations, A, and the pertur- 
bation parameter, A. 

We now substitute the decomposition of X ( 2.27 ) into the X-expansion of F ( 2.26| ) 
keeping only terms which will end-up being up to fourth order in A. Noting that X receives 
its first contribution at the second order^, and taking account of the parity symmetry for 
A (which will be justified for our case shortly) we find 



F(A, X; ii) = Fo(/i) + ^ (5/x A2 + F2{X, X) + A^ G{X) + F^gl A^ 



(2.29) 



"More explicitly, Fi{X,...,X) = J drdzFi{r,z) X{r,z) 
carrying derivatives. 



X{r,z) where here each X could also be 



^More concretely A is defined by A :=< X,Xgl > the projection of a general perturbation X onto 
Xql where < ■, • > is the appropriate inner product in the space of perturbations. As will become clearer 
in section M the inner product is given by the coefficient of k^ in the action, and makes the eigenvalue 
equation for the GL mode (|3.6|) self-adjoint. In the particular gauge used in this paper it is given by 



< X,X >:— const J aar"^ dr (3.21). The definition of X may be further described as follows. The 



spectrum of (3.6) has a single negative mode — fc^L ^nd s- continuum of scattering states X^ir) labelled 
by K := —k^ > 0. We may decompose the general perturbation X{r) according to the orthonormal basis 
of eigen-functions as follows X{r) = \XGL{r) + /.°° dhiX^n) X^ij-) where X{k) axe the coefficients in this 
basis. 

®5/i will turn out to be 0{\^) as well, see (^.12[). 
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This is the expression with which the LG analysis is carried out, and we now proceed to 
describe its various ingredients. ^ is a constant that can be read from dF/dfi. F2 is the 



same bi-hnear functional as in (2.26), F4GL is a constant given by 



Fagl ■= F/i{XGL, Xgl, Xgl, Xgl) 
= (O(A^)ofF(AXGL)) , 



(2.30) 



where by 0{\'^) we mean the coefficient of 0{X'^). Finally, G is a linear functional given 
by« 



G{X) ■.= 3F3{Xgl,Xgl,X) = 
'o(A4)ofF3(X,X,X)| 



X=AXgl+A2 



.)• 



(2.31) 



Substituting the perturbations series ( p. 23 ), or more explicitly ( 2.28| ), into the expan- 
sion of the free energy ( p. 29) ) we may solve for the back-reaction Xbr by varying the free 
energy. Equivalently, in practice we expand the Einstein equations to second order and 
obtain a source quadratic in Xgl- The equations of motion for the back-reaction, may 
be symbolically written^ as = 6F/SX = 2 F2XBR + G and its solution^*^ can be written 



symbolically as Xbr = —(1/2)^2^ G. 

Substituting back the fields into the free energy ( 2.29| ) one obtains the quartic coeffi- 
cient C as defined above in ( p. 10 ) to be 



C = F^GL - F2[BR) 



(2.32) 



where F/^gl was defined in ( |2.30| ) and F2{BR) is the quadratic action for the back-reaction 
defined by 



F2{BR) ■.= F2{Xbr,Xbr) = 

= (0(A4)ofF(A2XBR)) 



(2.33) 



Note that the negative sign in front of F2{BR) comes from the equation G = —2 F2 Xrr- 
Here we would like to note another more intuitive but less general interpretation of 
the criterion for transition order ( 2.12 , 2l3^ ). In cases where the original phase, before 
criticality, was stable, namely it was a local minimum of the free energy, the criterion for 
a second order transition is that the critical solution itself is also a minimum of the action 
(despite the presence of the zero mode), since clearly if it is not a minimum then another 
lower minimum already exists for a first order transition to occur. In such a case we may 
view the equation for the back-reaction as an attempt to find the direction of "steepest 
decent" and C would be the 0{\'^) coefficient of that decent. However, since we work in 



^Actually, since Xgl is a zero mode F2{X,X) = F2{X,X) for all X. 

^G(X) can be described more concretely in terms of the source for the back-reaction equations Src^ 
through G{X) = j dr Src^{r)X{r). 

^In vector notation it would read — 2F2,ijXgji + d. 

^"These back-reaction equations are solvable, namely F2{X,X) is invertible, since the zero mode Xql 
was removed from its domain of definition. 
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the canonical ensemble for which the (fat) uniform string has a negative Gross-Perry- Yaffe 
(GPY) mode [12|, the gravitational action is unbounded from below and the considerations 
above do not apply (or apply "modulo this negative mode"). 

Incorporating invariance under z -translations. The black brane is invariant under the U{\Y 
isometry group, originating from torus translations. As a result we can Fourier decompose 
all fields X. We can account for this decomposition in the computation of C thereby making 
( |2.32 ) more explicit. Here we discuss the case p = 1 (higher p will be discussed in section 



The U{1) isometry allows us to perform a Fourier decomposition of all fields as follows 

X{r, z)=Y, Xn{r) exp{i n kaLz) , (2.34) 

n 

where the index n is the harmonic, and in particular Xq is the z-independent mode, Xi is 
the GL mode etc.^^ 

The coefficient of the first order mode. A, is properly considered to be complex, as 
follows 

A X(i) = A e^ ^Gi ^ Xgpy (2.35) 

where we used Xql = ^^piikcL z) Xgpy |13|- It is seen that the phase of A yields z- 
translations. By translation symmetry all possible phases of A are equivalent. A natural 
way to fix the phase is to require A to be real. Choosing real A represents the spontaneous 
breaking of z-translations by the mode, but A — > —A is a residual symmetry corresponding 
to translation by a half period and this is the origin of the parity symmetry for A which is 
crucial for the expansion of the free energy as discussed above. 

The free energy is real, therefore U{1) invariance implies the following simplification 

F2{X,X) -^F2{Xo,Xo)+F2{X2,X2) + ... 

P G{X) -^ |Ap G{Xo) + [a2 G(X_2) + c.c^ , (2.36) 

where the indices of X denotes the harmonic, and the first equation is a result of the 
orthogonality with respect to F2 between different harmonics. Altogether the resulting 
free energy is 

F{X,X-f,) = Fo{fi)+A6fi\X\^ + F2{Xo,Xo) + F2{X2,X2) + 

+ |A|2 G{Xo) + \X^ G(X_2) + c.c.l + F^GL |A|^ , (2.37) 



where "c.c." stands for "complex conjugate". This form adds detail to ( ^.29 ). 



Since first order GL mode is (by definition) in the first harmonic, its square which 
sources the back-reaction has harmonics and 2, and therefore the back-reaction decom- 
poses into and 2 harmonics. Symbolically 

BR= {BRo,BR2) . (2.38) 



^^In the actual calculation we shall use a slightly different normalization, expanding into cosines and sines 
rather than exponentials, in order to account for X being real. 
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Finally 

C = F^GL - F2{BR^) - F2{BR2) , (2.39) 

is a more detailed version of ( p.32 ), where we made use of the orthogonality of different 



harmonics. 

3. A single compact dimension, T^ 

In this section we consider backgrounds with single compact dimension, p = 1 and D = 
d + 1. The most general metric in this case can be written as 

ds^ = e'^^f{r)dt^ + e'^^f{r)-^dr^ + 2Kdrdz + e^^ dz^ + e^^r'^d^l_2, (3.1) 

where the functions A^ B, K, H and C depend only on r and z. If they vanish and /(r) = 
1 — (ro/r)°'~^ then the uniform black string solution is reproduced, where tq designates the 
horizon location. (Below we set rg = 1). 

We are interested in constructing static perturbations about the uniform black string, 
hence the above metric functions will be considered as small corrections to the background 
metric. 

We must eliminate first the unphysical degrees of freedom by fixing a gauge. In our 
case, there are two degrees of freedom (related to diffeomorphisms of the (r, z) plane) that 
can be used to eliminate two of the five metric functions in (|3.1| ). It is not clear what is the 
"optimal" gauge choice and naturally we attempt here to simplify the equations as much 
as possible. We choose the gauge partially by requiring K = Q. The motivations and fixing 
of the remaining freedom is described shortly. 

3.1 First order — marginally tachyonic mode 

The small parameter of the perturbation theory. A, is the amplitude of the negative mode. 
To simplify the derivation of the equations of motions and expansion of the free energy in 
this section we use real A. This together with the U{\) symmetry along z suggests that 
the linear order perturbations are of the form 

X(^)(r,z) = Axi(r)cos(A:z), (3.2) 

for X = A, B,C and H (see also ( |2.35| ) ). After plugging these expressions into the Einstein 
equations, R^i, = 0, we obtain a set of ordinary linear differential equations (ODEs) for 
ai(r), 6i(r),ci(r) and hi{r). These equations will determine the first order perturbations 
and the critical wavelength kcL- 

At this stage we fix the gauge and our objective is to get the simplest equations 
(desirably uncoupled and without si ngular points between the horizon and infinity.) ^■^ 

^^One however should not be too overdihgent here: any gauge can be taken at the hnear order but at 
higher orders of the perturbation theory the equations might become degenerate, indicating that the gauge 
is too restrictive. For instance in our case one could choose H = 0. Then a GPY-sort equation mi h4, H] is 
reproduced at the linear order; it is easily solved in various dimensions, yielding kcL- Continuing with this 
gauge to higher orders of the perturbations theory leads to a contradiction as the back-reaction equations 
appear to put constraints on the first order perturbations. Physically, this is not surprising at all. Taking 
constant H is simply inconsistent with allowing non-uniform solutions, since when the non-uniformity 
develops the scalar charge defined by H must vary |lq, |l9[. 
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Examining the Rrz = constraint, 

2{d - 2)/ci + r/'ai - (2(d - 2)/ + rf) h + r {a'^ + {d- 2)4) = 0, (3.3) 

we are led to choose 

2{d-2)fC + rf'A 

^ - 2(d-2)f + rf' ■ ^'-'^ 



In particular, this is our choice at the linear order. From ( |3.3| ) we get then a[ = —{d — 2) c'^ 
and hence 

ci = -ai/(d-2), (3.5) 

as an integration constant, which could in principle arise, vanishes by boundary conditions 
at infinity. Consequently, the linear order equations are 



„ , {d-2)f + rf' 
J ai^ Oi + 



2 2{d-l){d-3)f 



12 1 



-k' + 



ai = 0, (3.6) 



[2{d-2) f + rff_ 



^d-2 y^ /^i; ^ 2{d-2) f + rf 



Besides, there is the constraint 



2/-rr , 2{d-l){d-3)f 

1 2{d-2)f + rf' 1 [2(d-2)/ + r/f ^ ^ 



We do not solve it but we check that it is indeed satisfied by the solution of (p^JST 

The equations are subject to boundary conditions: regularity at the horizon, that gives 

a' k"^ k"^ 

2{d-2), K = - ai at r = tq, (3.9) 



ai d — 3 d — 3 

and regularity at infinity, r ^ oo, which eliminates any growing solutions. 

We begin by solving (^]^). Because of linearity and homogeneity of (|3.6D the value 
ai(l) can be chosen freely. A particular value of ai(l) defines the normalization of the 
negative mode. We adopt 

ai(l) = -((i-2). (3.10) 

Then (3.5) gives ci(l) = 1, the normalization used in earlier works [0, H, 0]. 



Having fixed ai(l) we are left with a one-parameter shooting problem: only for a 
particular value of the parameter kcL ~ the critical GL wavelength - an integration outward 
the horizon converges. Adjustment of k and integration are iterated until k is found with 
desired accuracy. Our current implementation reproduces exactly the /ccL-values cited in 
the literature, see table |^. 

Next we solve for hi by integrating ( p.7p from the horizon outwards. A finite solution 
that asymptotes to a constant at large r's exists for any choice of /ii(l). Picking some 
/ii(l) we integrate, find the asymptotic value of hi, take minus this value for /ii(l) and 
reintegrate the equation to have hi vanishing at infinity. Namely, the shooting procedure 
for hi reduces to what we dub "second shot hits" . Clearly, the described procedure succeeds 
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because (3/7) contains only derivatives of hi so shifting the solution by a constant is still 
a solution of ( |3.7| ). This completes the first order computation. 

In closing we wish to stress the advantage of the gauge choice ( |3.4[) by comparing our 
master equation ( |3.6D for the negative mode with those considered in the literature. For 
example, in Gubser's gauge H, ^, which eliminates bi from (^^), two coupled equations 
for ai and ci must be solved simultaneously in order to find kc (the equations in an arbi- 
trary dimension are listed in Appendix of |^). Alternatively, in the GPY approach where 
one considers perturbations of the d-dimensional Schwarzschild solution in the transverse- 
traceless gauge, (which is exactly equivalent to setting H = 0) one obtains a single master 
equation for the negative mode [12, 14, 15, 0]. The equation, however, has an additional 



singular point apart from the horizon and infinity, whose presence complicates the shoot- 
ing^^(see Q for an analytical solution of this equation in a large d limit.) Recently, the 
Harmark-Obers gauge Q was considered in Q which obtained a single equation free of 
additional singular points. This is similar to what we got here. 

Having a single regular (aside from the boundaries) equation ( |3.6| ) we may arrive at 
the canonical form —d?ip/dp'^ + Vip = —k'^ip by changing he coordinate to p, related to r by 
dp/dr = /^^'^, and redefining ai to get rid of the first derivative. The equation mimics a 
Schrodinger-like problem of a particle with energy —k'^ moving in the influence of potential 
V. Figure Q depicts the potential and the wave- function of the negative mode in the case 
d = 5. 




Figure 4: The potential (solid line) and the wave-function (dashed line) for the negative mode, 
obeying the equation —d^ip/dp^ + Vtlj = —k'^ip, for d = 5. 



3.2 Back-reaction 

The back-reaction is comprised of the zero modes and of the second harmonic modes which 
we denote as 

X(2) (r, z) = A^ xo{r) + A^ X2{r) cos(2 k z), (3.11) 



^^ Basically, doing series expansion and matching from both sides about this point in a manner described 
in |15| allows solving the equation quite effectively. We have done this, verifying that the eigen-values kcL 



found applying this method are identical with those in table hi. 
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where X = A,B,C and H (see also (p^) 



As usual, in the presence of compact space the modes having dependence along the 
internal dimensions decay exponentially with r and they are massive in the Kaluza-Klein 
sense. The modes without z-dependence decay as inverse powers of r and they are massless. 

At any given order the gauge is given by (3.4). However, the gauge choice in higher 



orders can be modified by adding terms from the lower orders with the aim of simplifying 
the source terms. 

Massless modes 

Our method to simplify the equations is to choose the gauge such that the constraint 
Rrz = allows algebraic elimination of one of the fields. For the massless modes this 
constraint is trivially satisfied since these modes do not depend on z (actually Rrz is 
proportional to k times the LHS of ( |3.3D after replacing the subscripts 1^0). In this case 



for simplicity we adopt the same gauge as in (|3.4D 



2{d-2)fco + rf'ao 
'°= 2{d-2)f + rf' ■ ^'-''^ 

Note that the variation of temperature at the horizon, which is proportional to exp(A — B), 
vanishes automatically in this gauge. 

The equations governing the massless mode can be written schematically as 

^-id-2) (^r^-^fh'^"^' = Src{a,,a[,h,,h[) 

£a{aQ,aQ,c'Q,ao,co) = Src{ai,a[,hi,h[;ho,h'Q) 

£c{co,aQ,CQ,ao,co) = Src{ai, a[, hi, h[; ho, h'^) 

Q(ao,Co,ao,co) = Src{ai,a[,hi,h[;ho,hQ), (3.13) 

where the sources contain squared first order perturbations, the exact form of the operators 
£a,£c,Q and of the sources is found in Appendix ^. The last equation is a constraint - it 
is not solved but we verify that it is satisfied by the solution of other equations in ( 3.13| ). 



The equations ( 3.13 ) are subject to regularity boundary conditions at the horizon 
( |A.5| ). These conditions determine the derivatives of the functions in terms of ao(l), co(l) 
and /io(l) which are 3 free parameters. At infinity we demand the length of the compact 
circle and the period of the Euclidean time to reach their unperturbed values, namely we 
demand that A = H = Oasr^oo. Besides, an obvious requirement is that C approaches 
a constant at large r. However, it turns out that any choice of that constant yields regular 
solution to ( |3.13|) . In addition, the constraint is satisfied for any asymptotic value of C 
and it does not contain any additional information beyond that already known from the 
second order equations. Therefore, we end up with 5 out of 6 conditions/parameters which 
are necessary to specify a unique solution to the three second order ODEs in (|3.13| ). 

The situation is puzzling only until one realizes that the choice ( 3.12| ) does not fix 



the gauge completely. By analyzing the residual gauge we find that fixing the residual 
reparameterization of r is equivalent to choosing the asymptotic of cq. To examine the 
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effect consider a shift 

r^r + ^{r), (3.14) 

that induces the variation in the metric functions: 

^ao = i^^, Sco = -, 5ho = -l-A + i', 6ho = 0. (3.15) 

2/ r 2/ 



Substituting this into ( 3.12| ), and demanding invariance of the gauge one gets a first order 



ODE for ^. The equation has a solution such that ^(1) = and that asymptotically 
C = Coo'^; with some constant ^^o- It follows from ( 3.15| ) that ^oo is the asymptotic value 



of cq. Altogether we fix this remaining freedom by requiring C = at infinity. ^^ 

To solve ( |3.15 ) we first treat /iq, which is decoupled from other metric functions. 



The solution is achieved in the "second shot hits" fashion since the equation contains 
only derivatives of ho, similarly to what occurred at the first order. Then we solve two 
coupled equations for oq and cq where ho as well as squared first order perturbations act 
as sources. This is a two-parameter shooting problem: the values co(l) and ao(l) are 
adjusted and the equations are integrated from the horizon outwards iteratively until a 
solution with asymptotically vanishing ao and cq is found. This completes the massless- 
modes computation. 

Massive modes 

The equations for these modes are similar to those we had at the linear order (3.6-3^) 



but this time there are sources. To see if we need to modify the gauge (3^) by adding 
some lower order terms squared we examine again the Rrz = equation: 

2{d - 2) f C2 + r f a2 + 2 {a'2 + {d - 2)4) + 

^W^'' ^ "" "' - (^ (^ - 2) ^ + ''^') ^2 = 0- (3-16) 

From here the obvious choice is 

_ 2{d-2)fc2 + rfa2 d-l rfa,a[ 

' 2{d-2)f + rf ^2(d-2)2((i-2)/ + r/'' ^^ ' 

since it allows algebraically to eliminate one of the fields, in precisely the same manner as 
at the linear order. Consequently, C2 = —a2/{d — 2) and so we are left with only two second 
order ODEs to be solved. 

Schematically, the equations read 

8a{a2,a'2,a2) = Src{a[,ai,h[,hi), 

£h{h2,h'2) = Src{a'i,ai,h[,hi;a2) (3.18) 

and their exact form is listed in ( [A.?] ). These equations are subject to regularity boundary 



conditions: both at the horizon, ( A.9| ), and at infinity. 

'^^A similar phenomenon was encountered by Gubser M who fixed the residual gauge, remaining after 
making the conformal ansatz, by assigning co(l) some arbitrary value. Different choices of this value set 
different "schemes" in Gubser's language. 
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The 02 equation is independent of /12 and it is solved first. To this end we shoot with 
02(1) as the shooting parameter. It is adjusted until the regular (decaying at large r) 



solution is found. The equation for /i2 contains 02 as a source, see ( 3.18| ) or (A/7), and 
once 02 is known /i2 is found by the "second shot hits" method. Finally, we verify that the 
constraint ( [A. 10 ) is indeed satisfied by the solution. 



3.3 Free energy and comparison ^vith previous results 

As soon as the negative mode and the resulting back-reaction are found we are ready to 
compute the coefficients of the free energy expansion about the critical point up to the 
fourth order in A. 

The quadratic term in the expansion ( |2.29D , F2 {X, X) , is obtained by expanding the 
York-Gibbons-Hawking llB, 17| action integral (2.2C|) in our gauge. We find 



j(2) 



IQttGd 



dz 



drr 



d--2 



ro 



f{dra{r,z)f + {dMr,z)f 



d-1 
2{d 



l)(d-3)/ 



/2 



■ a{r, zY 



. ■. . , (3-19) 

{2{d-2)f + rfr 

Note that h{r, z) does not appear here. To get J''^' in terms of harmonics an{r), such 
as ( p!^ , |2.38D one needs only to substitute in the harmonic decomposition^^ a{r, z) = 
Yln=o^n{r) cos{nkz). 

The quadratic part of the action ( 3.19| ) allows us to compute A ( |2.29 ). By definition 
the integrand vanishes for k = kci and ai(r) being GL mode (found in subsection ^H] ); 
the leading term appearing as one moves away from kcL, reads 



/(2) 



P^d^2 

16-kGd 



dz 



ro 



2kGLSk 



d-1 

2{d-2)' 



„d-2 



aiirf 



dry 



(3.20) 



Using the relations (|2.1| , |2.7D and (^ 
as 5k = 2/{d — 3)6iJ,. It follows then 



we connect the dimensionless wavenumber with /u 



A^F^^^/X^ 



oc 



47r(d-l) 



a? r'^-^dr 



(3.21) 



{d-2){d-3) ji 

The factor of proportionality, Qd-2/{lQ'^GD), appears in all action integrals such as ( 3.20| ). 
So in concrete calculations we shall get rid of it simply by suitably redefining the coefficients 
of the free energy expansion. 

The quartic contribution to F is ( ^.391) . We omit the explicit expressions for the inte- 
grand of F/^GL since it is far more cumbersome then ( 3.ig| . 3.21 ), though it is straightforward 
to obtain. We note that in our gauge (|3.4|) F4GL includes pieces which originated in the 
metric function B before the gauge was fixed. ^^ 

The numerical values of A and of the specific terms forming C are listed in table |2|. It 
is evident from this table that a critical dimension, D*^^ = 12.5, for the canonical ensemble 
appears, above which C changes sign marking a change in the phase transition's order. 

^^ Actually, we checked this formula only for n 7^ 0. The zero mode sector may contain additional terms. 

^''It is also possible to carry the expansion of the action before gauge fixing, which was the way our 
computation actually happened to be performed and only later was it translated into a decomposition of 
the gauge-fixed action. The bottom-line result for C is of course the same in both ways. 
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d 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


A 


0.465 


2.05 


5.11 


9.98 


16.9 


26.2 


37.5 


51.7 


68.7 


88.6 


112 


C 


-0.403 


-1.39 


-2.99 


-5.03 


-7.07 


-8.47 


-8.30 


-5.27 


1.97 


15.0 


35.7 


Fagl 


.0976 


.7883 


3.035 


8.417 


19.27 


38.85 


71.42 


122.48 


198.8 


308.8 


462.2 


F2{BRo) 


.3918 


1.287 


2.239 


1.932 


-2.112 


-13.82 


-38.86 


-85.06 


-162.4 


-283.4 


-463.3 


F2{BR2) 


.1076 


.8912 


3.791 


11.52 


28.46 


61.14 


118.6 


212.8 


359.2 


577.2 


889.9 



Table 2: Coefficients of the free energy expansions. We notice that C changes sign between d= 11 
and 12. This indicates that in the canonical ensemble a second order phase transition occurs in 
dimensions d> 12. 



Landau-Ginzburg method, derived from A, C 



d 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


0-2 


1.40 
-0.569 


3.00 
-0.976 


4.68 
-1.40 


6.21 
-1.79 


7.36 

-2.08 


7.93 
-2.21 


7.68 
-2.09 


6.29 
-1.70 


3.56 
-0.969 


-0.695 
0.197 


-6.76 

1.85 


Previous method, direct computation 


d 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


0-2 


1.45 
-0.595 


3.03 
-0.989 


4.69 
-1.42 


6.21 
-1.80 


7.34 
-2.02 


7.89 
-2.14 


7.66 
-2.08 


6.25 
-1.66 


3.57 
-0.966 


-0.724 
0.198 


-6.84 
1.82 



Table 3: Coefficients of the dimensionless mass and of the entropy variation as defined in ( 2.18 2.19| ) 
computed in both the LG method and the previous method. The match is within 5%. The sign 
change (in both quantities) between d — 12 and d ~ 13 indicates a change in the order of the 
transition in the micro-canonical ensemble. 



There is a subtlety involved in the calculation of C. According to ( ^.32 ) it is defined 
as a difference between two numbers, -F4GL and F2{BR). As the dimension increases 
these numbers grow and they become comparable around D*^^, see table ^. Hence, it is 
imperative to compute each term in ( 2.32| ) accurately enough in order to obtain reliable C. 



Having computed the coefficients in the free energy expansion we may obtain the micro- 



canonical thermodynamic coefficients rji and a2- These are defined in ( 2.142.18 , 2.19 ) and 
their numerical values are summarized in table y. We include for comparison the same 
quantities computed in [^ (see also @) in Gubser's method. We do not know which 
method is more accurate, but obviously both yield comparable results (with less than 5% 
discrepancy.) 

Notice that for each dimension both methods produce exactly two numbers that encode 
the entire thermodynamics at leading order: 771 and C72 in the micro-canonical ensemble 
or A and C in the canonical ensemble, and we have demonstrated the relation between 
the two pairs of quantities. We note also that while in the canonical ensemble ui = 
automatically, in the micro-canonical ensemble this is a derived result. In fact, in that case 
the smallness of ai can be used as an estimator of numerical error B. 
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4. Several compact dimensions 

After having demonstrated the LG-inspired method, we would hke to apply it to the 
computation of the order of the transition in the presence of a more general compactification 
- by a torus. 

Square torii. We choose to restrict our study to square torii, namely z* ~ z* + L with 
the same L for all z'^ and with right angles between the axes. We motivate this restriction 
by the following. We view the space of torii as having two boundaries - on the one hand 
highly asymmetrical torii, where one (or more) dimensions are much larger than the rest, 
and on the other hand highly symmetrical torii such as the square torus. 

For example, for p = 2 the space of torii is the well known modular domain parame- 
terized by r. The highly asymmetrical 2-torus boundary is to be found at r — > ioo. The 
most symmetrical torii, the square and hexagonal ones are to be found on the boundary at 
r = i,exp(i7r/3) respectively. The boundary of the modular domain also includes other, 
less symmetrical torii. Actually, we could have chosen to study the hexagonal torus - it 
also enjoys a large symmetry that would simplify the calculations, and our choice is simply 
one of convenience. 

The limit of highly asymmetrical torii is easy to understand. Suppose one of the edges 
of the torus is much longer than the rest (or equivalently, one of the vectors in the reciprocal 
torus is much shorter than others). In this case as the mass of the brane is reduced the 
first GL instability to occur will be associated with this long direction. Since this mode 
is invariant under translations in all other directions, and as long as that symmetry is 
not spontaneously broken (actually it may very well break to some extent during the 
collapse when the highest curvatures develop) then we find that we can completely ignore 
all directions except for the largest one. So effectively that case is reduced to that of a T^ 
compactification. A similar reduction applies if there are several directions which are much 
larger than the rest - in such a case we may reduce the problem to one with a smaller p. 

Therefore by studying the limit of the symmetrical, square torus, we hope to achieve 
an understanding of both limits and thereby also some understanding of the intermediate 
region of general torii. 

GL instability. In a T^ compactification (this discussion appeared already in [y) we have a 
GL instability for each vector in the reciprocal lattice. This replaces the single instability 
and its harmonics in the T^ case. Therefore as the mass of the black brane is reduced these 
instabilities appear in order of smallness of the vectors in the reciprocal lattice. Of course, 
once the smallest one is encountered the system will start collapsing and we will not have 
a chance to observe the other instabilities separately. 

For symmetrical torii, including the square, a non-generic phenomenon happens: as 
the first instability is reached (say by lowering the mass) several modes turn marginally 
tachyonic simultaneously. More precisely, for T^ there are p such modes. One may be 
concerned that this degeneracy occurs only for exactly square torii, but we would like to 
argue that this degeneracy is relevant also for torii which are nearly square. Indeed, if 
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the torus is not a precise square one of the GL instabihty modes will be triggered first. 
However, as the other "would be" tachyons have a very shallow potential at this moment, 
once the instability starts "rolling" and energy becomes available there is nothing to stop 
these modes from getting spontaneously excited, thereby allowing motion of the system in 
all the p "tachyonic" directions just as for an exactly square torus. 

Generalizing the perturbation method. We would now want to generalize the perturbation 
method and find a generalization for the coefficient C ( 2.32D which determines the order of 



the transition. The first order perturbation is an arbitrary linear combination of the p GL 
modes, with coefficients denoted by Aj, namely the generalization from p = 1 to arbitrary 
p is given by 

X(i) = XXgl^ X(^^ = Xi X}.^ , (4.1) 

where indices i^j are to be summed in the range 1 < i,j < p and Xqj^ is the GL mode 



associated with z* (more explicitly, the mode computed in subsection 3.1 after substituting 
z -^ z*). 

Proceeding to the second order the sources have the following schematic form Src ~ 
Xi Xj Xq^ ^gl ■ Accordingly the back-reaction becomes 

Xbr = Xi Xj X'^j^ . (4.2) 

To see the implications of the symmetries of the square on X^^j^, let us consider an 
arbitrary abstract tensor T*-'. Since the symmetries of the square allow us to exchange 
Zi ^^ Zj for any i and j it is clear that there are only two independent components to T*-' 
depending on whether i = j oi i ^ j, which we denote as T^ and T^ respectively. In 
summary 

T"=[i; i;i ^ (4.3) 

We shall also make use of another definition, "T-average" 

- _ pT=+p{p-l)T^ ^ T= + {p-l)T^ ^^^^^ 

Thus there are two kinds of back-reaction. BR^ is the one that was computed already 
fQ],i7 p = i^ while BR^ is novel and will be determined in the next section. 
Now we must substitute the perturbation series 

X = XiXQj^ + Xi Xj X^p^ (4.5) 

into the action and compute the term of 0{X^). Concentrating on the quartic contribution 
of the GL modes we find the following generalization 

F,gl\X\^^FI^^\X,\^\X,\^ (4.6) 



More precisely, for each i, BR" is the mode computed in subsection 3.2 after substituting z 
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(C/(l)^ invariance requires the terms to have equal amounts of A* and A*). Similarly, the 
other factor in the C formula ( |2.32| ), F2{BR) gets generalized as 

F2{BR) |A|^ ^ F2{BR'\BR") |Af |A^f + F2{BR'^ ,BR'^) |Af jA-'f (4.7) 

(note that due to orthogonality of harmonics only the 0-harmonic contributes to F2{BR^^, BR^^) 
with i 7^ j). 

Altogether we find that C is replaced by a tensor 



Determination of the order. Consider C(A*) as it varies over all possible directions in 
tachyon space A*. It is enough to consider unit vectors X] |Ajp = 1, and thus the vector 
I Aj p varies over the simplex given by this normalization condition together the inequalities 
< |Aip (and < 1). We shall now show that the range of C(A') is precisely 

C(A') 6 [C-,C(p)l c [C-,C*\ (4.9) 

V 

where [a, 6] denotes here the interval between a and h irrespective of which one of a, 6 is 
bigger. To see that we compute 

C(A') =C^\\\^\\jf = 

=c-j:iA.r+c^((5:iA.if-j:iA.i^) = 

= C^+(c=-C^)^|Ai|^ (4.11) 



Since ^ |Aj|^ ranges precisely over [1/p, 1] ([4.9D follows. 

The transition is second order if and only if C(A*) is positive for all directions in 
tachyon space A*, for otherwise, if there exists a vector A* such that C(A*) < the system 
will spontaneously settle on that direction and a first order transition will ensue. We 
conclude that the transition is second order if and only if both 

C=, C{p) > 0. (4.12) 

p-dependence from p = 2. From their definition it is evident that C^, C^ are p-independent. 
Actually C^ is already known from p = 1 

C==C (4.13) 

Moreover, it suffices to compute C^ for p = 2 where the mixed ij terms appear - the 
computation does not change for higher p. 
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Thus it remains to determine C^. Rather than determine it directly we take a some- 
what more physical approach and pick the diagonal direction in A space, namely 

A;=^A (4.14) 



independently of i. C{\i) computed along this direction is precisely C defined in ( 4.1C| ). 
Therefore by computing C|p=2 we obtain 

C^ = 2C\p=2-C= . (4.15) 

The A term in the free energy may be easily generalized as follows 

AdfilXl"^ ^A6fiX'~X\ (4.16) 

and on the diagonal direction ( [4.14] ) A = A. 

The actual computation (at p = 2), which is fully described in the next section, requires 
obtaining the mixed term F^q^ l^iP l-^2p from substitution into the action, a new mixed 
term in F2{BRq), computing the mixed back-reaction term BR"^ = BRn, which is in 
the (±1, ±1) harmonics, and finally substituting it into the quadratic action to obtain 
F2{BRii). All the other ingredients were essentially computed already for the p = 1 case. 

Micro- canonical ensemble. Having obtained the "diagonal" coefficients C and A we can 



compute the variation in dimensionless mass r/ ( 2.17 ) along this direction by using ( 2.18 ) 



nl■■=-'2{d-?>)-^ + -^^iGLA (4.17) 

^GL A d-2 



In fact, we can define a tensor of mass variation 771, as in (|4.4| ). Then, in analogy with what 
we had for the tensor C, we have r/f = ??i|p=i, and recalling ([4.15| ), rf[ = 2fji\p=2 — Vji ■ 
By analogy with ( [4.12 ) the transition will be second order in the canonical ensemble for a 



given p exactly when both 

r?r,f?i(p)<0, (4.18) 

where as before ( [4.10| ) we define r^i(p) := Vji /p + (p — 1) Vi /P- 

5. Computation for T^ 

In this section we perform the computation for the square two-torus, so D = d + p with 
p = 2. 

The most general metric consistent with the U{l)t x SO{d — l)n isometrics is 

ds^ = f{r)e^'^dt^ + f{r)-^e^^dr'^ + Kidz'dr + Hidz'dz' + Wdz^dz^ + r'^e^^dnl_2 (5.1) 

where i = 1,2. The metric in the (r, z) space has 3 gauge degrees of freedom. We fix 2 of 
them by choosing Ki = 0. The remaining degree of freedom will be fixed by considerations 
similar to the T^ case. 



23 



The square torus has also the fohowing three discrete isometries, which wih also be 
preserved by our perturbation 

Zi ^ Z2 

Zi -^ -Zl 

Z2 ^ -Z2 (5.2) 

First order — marginally tachyonic modes 

In the square torus background two GL modes turn marginally tachyonic simultaneously, 
resulting in a two-dimensional "tachyon space". In this section we focus on the diagonal 
direction in the "tachyon space" , given schematically by Aj X^j^ where 

Yi = \ (5.3) 

Here it is convenient for us to use this normalization to compute the back-reaction and 
turn to the normalization ( [4.14| ), which differs by a \/2 factor, later when we substitute 
into the free energy (which amounts to dividing the quadratic and the quartic coefficients 
of the free energy hy p = 2 and by p^ = 4, respectively). 

As in T^ case we use real A. The diagonal instability is given explicitly in terms of the 
functions ai(r), bi{r), ci(r), hi{r) and the critical constant k = kci which were found in 



subsection 3.1 as follows 



for X = A,B,C and 



(r,z) = 


Axi(r) y^( 


20s{k z 


Hi{r,z) 


= \hi{r) cos{k Zi), 


W 


= 0. 





(5.4) 



(5.5) 

Clearly, because of the symmetry between the tachyons along the torus edges only a single 
function hi{r) is used for both i's, and actually the relation zi <-> 2:2 ^ Hi <^ H2 will hold 
throughout the perturbation, so essentially the H^s reduce to a single function. Note that 
the discrete square symmetries (|5.2D are indeed preserved by the perturbation. 

As expected, after plugging the expansions into the Einstein equations, the equations 
for different z's completely decouple. Besides, the gauge conditions and the equations are 
identical for every i and they coincide with those derived in the T^ case, see ( p.6| - |3.8D . The 
upshot is that the first order computation in the T^ (T^ in fact) case is essentially identical 
with T^ 

Back-reaction 

The form of the back-reaction sources can be roughly obtained by squaring the linear order 
expansions (5^) and (|5.5| ). In addition, a mixed W term, which was absent in the T^ case. 



becomes possible. 
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Massless zero modes 

The expansion here is essentially the same as for T^. 

Xir,z) = 2A2xo(r) , 

W = , (5.6) 

where X = A,B,C and the p = 2 factor comes from summing over the back-reaction 



zero-modes for all i. The vanishing of W is dictated by the square symmetries (5.2). 



To confirm this form we proceed as follows. We use the gauge ( ^.12 ). The ODEs 



obtained after substituting the above expansion into the Einstein equations are given by 



( A.1 , A.2 ). The equations differ from the T^ case only by the appearance of the factor 
p = 2 that multiplies the sources in ( [A.2| ). As a result, ho{r) is unchanged while ao and cq 
are modified - multiplied by p = 2 relative to T^ solutions. Indeed, that was the reason 



that we introduced the factor of 2 in (5.6). 



Massive second harmonic modes 

The treatment of these modes replicates the one we had for p = 1. The obvious 
expansions 



X{r,z)=\^X2{r)^cos{2kzi), (5.7) 



i=l 



for X = A,B,C and 



Hi{r,z) = A /i2(r) cos(2 A; Zj), 

W =0, (5.8) 



yield exactly the same equations as for the one-torus ( |A.7| ). So nothing new needs to be 
solved. 

Massive mixed modes 

This is the longest section here, for these modes are the only new ones that were not 
discussed in the T^ section. 

The only mixed modes expansions allowed by the symmetries ( |5.2| ) are 

X{r,z) = A xii(r) cos(A;zi) cos(/!; Z2), (5-9) 

for X = A,B,C,H and 

W{r,z) = X wii{r) sm{k zi)sm{k Z2)- (5.10) 

Plugging this into the Einstein equations we examine the constraints R^Zi = which take 
the form 

2{d-2)cuf + rauf' 2 (d - 1) rfaia[ 



m 



2{d-2)f + rf' d-2 [2{d-2)f + rf'\ 

2rf {2f-rf')aih[ f r [2h'n + w[^ + 2a[^ + 2{d - 2) c'^] 
[2{d-2)f + rff 2{d-2)f + rf' 



0. (5.11) 
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As usual, this equation instructs us how to choose the gauge and we take 
^ 2{d-2)cuf + raur 
'' 2{d-2)f + rf> 

2(d-l)r/aia; 2r f {2 f - r f) a,h[ 

-^ {d-2)[2{d-2)f + rf']^ [2(d-2)/ + r/f ■ ^' ' 



This gauge is consistent with our basic choice ( p.4D which at this order only gets modified 
by adding sources from lower orders. 

Now ( ^.lll) reduces to 2 h[^ + w[i + 2 a[i + 2{d — 2) d^^ — which we solve for h\\ 

/ill = - 2 Ki + 2 an + 2 (d - 2) cn] , (5.13) 

where the asymptotic boundary conditions are used to eliminate the integration constant. 
The rest of the equations are given by ( A.ll|) . Their schematic form is 



'?a(oii,aii,Cn,aii,cii) = 5rc(ai,ai,/ii,/ii) 
'^c(cii,Cii,aii,cii,aii) = 5rc(ai,ai,/ii,/ii) 

r-^'^'') (r'^-'Mi)' = 5rc(ai,a'i,/i'i;aii,cii), (5.14) 

These equations are only partially coupled and can be further separated. The first two do 
not involve w\\ and moreover one gets a decoupled equation for a\\ — c\\. Having solved it, 
we proceed to solve first for a\\ + cn and then for w\\ where at each step earlier solutions 
appear as sources. 

The equations are subject to the horizon boundary conditions ( |A.13 ) and asymptot- 



ically they must vanish. We shoot to solve the equations, and the values of the shooting 
parameters are given in table 0. Finally, the constraint ( A.14| ) is verified to be satisfied. 



Free energy and thermodynamics 

In the rest of this section we change the normalization from ( |5.3| ) to ( [4.14 ) 



The quadratic term in the free energy expansion remains unchanged relative to the T^ 
case: Ap=2 = A = Ap=\. 

We turn to the quartic term 

C = F4GL-[i^2(5i?0)+i^2(5i?2)+i^2(5i?ll)] • (5.15) 

We can figure out which additional terms appear in the zero mode contribution for T^ 
relative to T^ 

F2{BRo)\p=2 = F2{BRo)\p=i --I r'^'^ j^>^2^^_ ^^_^^^ 

The second harmonic contribution requires only a p-dependent factor (division by p = 2) 

F2{BR2)\p=2 = ^F2{BR2)\p=i. (5.17) 

The mixed modes term, F2{BRii), is new, and likewise there is a new contribution to 
FiGL- It is straightforward to substitute the expansions ( ^.9| , ^.10 ) into the action integral 



(remembering to divide by 4 due to change of normalization). The resulting expressions 
are cumbersome and not particularly illuminating; we omit their explicit form. 



The numerical values of the various terms in (5.15) are recorded in table 
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d 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


c 


-.441 


-1.27 


-2.28 


-3.24 


-3.76 


-3.30 


-1.23 


3.17 


10.9 


23.2 


40.6 


Fagl 


.2869 


2.403 


9.864 


28.97 


69.59 


145.9 


277.5 


489.4 


813.9 


1290 


1967 


F2iBRo) 


.2775 


.8456 


1.130 


-0.2807 


-5.992 


-19.99 


-48.14 


-98.35 


-180.7 


-307.8 


- 495.4 


F2{BRii) 


.4023 


2.377 


9.121 


26.76 


65.24 


30.57 


59.29 


478.2 


804.0 


1287 


1977 



Table 4: The non-trivial coefRcients of the free energy expansions for T"^ . The rest of the coefficients 



are found using the fact that A is unchanged relative to T^ case, equation ( 5.17 ) and table 



d 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


Cp=2 


-.441 


-1.27 


-2.28 


-3.24 


-3.76 


-3.30 


-1.23 


3.17 


10.9 


23.2 


40.6 


c= 


-0.403 


-1.39 


-2.99 


-5.03 


-7.07 


-8.47 


-8.30 


-5.27 


1.97 


15.0 


35.7 


c+ 


-.481 


-1.14 


-1.57 


-1.44 


-.453 


1.87 


5.83 


11.6 


19.9 


31.3 


45.5 



Table 5: Components of the tensor C (|4.8| ), that determines the order of the phase transition in 
the canonical ensemble. 



d 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


m,p=2 


1.49 


2.81 


3.89 


4.65 


5.02 


4.84 


3.99 


2.40 


-0.125 


-3.73 


-8.42 


vT 


1.40 


3.00 


4.68 


6.21 


7.36 


7.93 


7.68 


6.29 


3.56 


-0.695 


-6.76 


vt 


1.58 


2.62 


3.09 


3.10 


2.68 


1.76 


0.299 


-1.48 


-3.81 


-6.77 


-10.1 



Table 6: Components of the tensor 771 (2.17), that defines the dimensionless mass variation (see 
definition at the end of section 0) which determines the order of the phase transition in the micro- 
canonical ensemble. 



p-dependence 

Table g summarizes our results for the C tensor. For the various values of d, the number 
of extended dimensions, we give C = C^ which is the p = 1 value, together with C\p=2- 
From these we determine C^ using ( [4.15 ). C^ allows to extend the results to any p by 
computing C{p) from ( [I.IOD , which in turn determines the order of the transition in the 
canonical ensemble through ( 4.12| ). 

For the micro-canonical ensemble C is replaced by the r/i tensor, and table ^ summarizes 
the results, which were obtained through the use of (2.18,4.17). The condition for a second 



order transition is given in (4.1 



6. Discussion of implications 



We demonstrated that the Landau-Ginzburg inspired method to compute the order is some- 
what simpler than the previous one, and gives the same results (tables ^^. In addition 
we analyzed the order of the Gregory-Laflamme transition for (square) torii compactifica- 
tions. The main results are summarized in tables |^|^ and we now proceed to discuss their 
implications. 
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Transition order is independent of p. A necessary condition for a second order phase 



transition is < C^ = C (4-12). Therefore, for all D < D*^^ = 12.5, where the transition 
is first order for p = 1 (namely C < 0) it is first order for all p. 

From table |^ we find that the converse it true as well. For d > D*^^ — 1 = 11.5 where 
C > then also C^ > and therefore from ( [4.10D also C{p) > for all p, and hence the 
transition is second order in this range of d for all p. 

Altogether we conclude that the transition order depends only on d, the number of 
extended dimensions, and not on p, the dimension of T^. It would be interesting to know 
whether this property holds for a general compactification manifold, Y^. 

Failure of the "equal- entropy for the equal-mass estimators" . Our original motivation for 
exploring the phase transition order in torus compactification arose after observing in Q 
that the "equal entropy for equal mass" estimator indicates lower critical dimension for 
some higher-dimensional torii (lower then the p = \ value of D* = "13.5"). In fact, this 
argument gave an encouraging value of "10" for T^ with 3 < p < 6. However, as just noted 
the critical dimension cannot be lowered, and here we explain why the estimator failed. 

We shall now show that the estimator is strongly sensitive to the approximation error 
which is inherently built into it. In fact an approximation error in the percent and sub- 
percent level will even cancel the prediction of a critical dimension altogether. In other 
words, it turns out that the estimates of H carry unusually large error bars, which resolve 
the puzzle. 

The estimator is constructed by comparing entropies of a uniform black-brane and of 
a localized black-hole with the same mass. The idea is that if entropy equality is achieved 
below the critical mass this may indicate that the phase transition between these objects is 
second order since a new non-uniform branch emerging from the GL point is likely to exist 
and to connect between the branches. The dimension (not necessarily integer) for which the 
"equal-entropy mass" coincides with the instability mass estimates the critical dimension. 
This is only an estimate because the entropy of a localized solution is approximated in this 
approach by a naive Schwarzschild value. This reasoning worked remarkably well in the 
T^ compactification ||3[, producing a fairly precise estimate of D* . 

However, the Schwarzschild formula is only an approximation to the localized black 
hole entropy, hence we use Sbhiv) = Sschwi'n){^ + f)i to model the modified expression. 



Perturbation theory pQ , 21] and numerical results [^, ^] indicate that the entropy of a 
localized black hole is higher than that of a Schwarzschild black hole and hence we must 
take e > 0. 



Solving for r/ we get 
1 



'75(e) 



16 vr 






{l + ^){D-Z)(d--i)/p (g^i) 



To find the estimate for the critical dimension, D* we compare this with the critical mass, 
which is well approximated by the formula |y], r]c{D,p) = r]c{D — p) = 0.47 • O.GSG^"^"*"^. 
Exploring the resulting dependence of this estimate on e shows that D* increases fast as it 
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grows. For example, in a p = 1 case, while for e = one has D* ~ 12.5, taking e = 0.002 
results in D* ~ 16, and e > 0.003 does not have a critical dimension at all. The situation 
is similar for other p's and for p = 6, e ~ 0.08 sufficies^^ to send D* — > cxd. 

In conclusion, there is no contradiction between our present results and the estimate 
due to the sensitivity of the latter to the actual value of the entropy of the localized phase. 

Subtle spontaneous symmetry breaking. Examination of table |5| reveals that for all^^ d > 4 
C\p=2 > C therefore the diagonal direction in tachyon space is disfavored relative to turning 
on a "single tachyon" namely a tachyon which depends only on a single z*. We interpret 
that to mean that the discrete symmetries of the torus are spontaneously broken: for a first 
order decay the time evolution will proceed (at least initially) through a single tachyon, 
while for a second order transition the system will re-settle into a slightly non-uniform 
string where mostly a single tachyon is turned on. This makes sense if we recall that the 
square torus is a special torus, and for all nearly square torii the degeneracy of the GL 
modes is removed anyway (see section ^) . 

Micro- canonical ensemble. From table |^ we find that the behavior in the micro-canonical 
ensemble is similar to the case of the canonical ensemble (except for the known change in 
the critical dimension). Since rj^ < vj^ for all d > 4 in our table, then j]"^ alone determines 
the order of the phase transition, and the order is independent of p (for all d). Moreover, 
as discussed in the previous paragraph we conclude that a single tachyon is preferred over 
the diagonal direction through spontaneous symmetry breaking (for all d > 4). 
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A. Second order perturbation equations 

Here we give the full back-reaction equations for each harmonic, including the full expres- 
sions for the sources. 

Zero-modes 



^*We see that e grows as p increases. This can be attributed to the fact |24[ that for a given d the 
corrections to Schwarzschild entropy is larger for larger p's. We thank Dan Gorbonos for sharing his 
unpublished results regarding small mass corrections for caged black holes in T*" compactifications. 

^^The case d = 4 is marginal. The numerical values of C\p=2,C are close enough, such that if both of 
them had a 5% numerical error (see discussion in section H) the inequality would hold for d = 4 as well. 
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First we solve the equation for ho 



il {r'''fh'o 



k^ hi 2{d- 1) {d - 3) /' a[ 



+ 



2/ [2{d-2)f + rf>\ 



2f-rf 
2 /2(d-2)/ + r/'"^ 



4(d-l)2(d-3)2/' 2 



f/21 



fe2[4(d-2)/ [(d2-3d+l)/ + dr/']+rY'j 2 ,ai^ 

ai- (A.l) 



[2((i-2)/ + r/f ' 2{d-2)f[2{d-2) f + rff 

Then the equations for oq and cq are solved 

2{d-2ff + 4{d-2)rff' + r^ f'^ 



a'/.+ 



■an + 



° rf[2id-2)f + rf'] « 

(d-2)/' [2(d-3)/ + r/'] , , (d_2) (d-3)/'2 



■Co + 



c'o' + 



2/[2(d-2)/ + r/'] " ■ /[2(d-2)/ + r/'] 

2 (d - 2)(2 d - 5) /2 + 4 (d - 2) r / /' + r2 /'2 



r/[2(d-2)/ + r/'] 
2 /' (d - 2)(d - 3) /2 + (2 d - 5) r / f + r^ /'^ 



■Cn + 



2 (ao -Co) = 5rCao, 
2 (d - 2) 



g'q 



rf[2{d-2) f + rff 
where the sources are given by 

fc2 (r /' - 2 /) aj l^hl ^ 2{d- 1) (d - 3) /' ai' 



r[2{d-2) f + rf] 
{ao - Co) = Srcco, (A.2) 



Src, 



ao 



2f[2{d-2) f + rf] 



2/ 



+ 



Src. 



CO 



fc2/ii 2(d-l)(d-3)/'ai' 



+ 



2(d-2)/ (d-2) [2(d-2)/ + r/'] 



[2(d-2)/ + r/r 



Ol 



2/ 



,(A.3) 



^ + 



2(d-l)2(d-3)f2 + (d-2)(2/-r/0[2(d-2)/ + rf]fc2 ^2 
2(d-2)2/[2(d-2)/ + r/f ""' 

The above equations are subject to regularity boundary conditions at the horizon 

fe^ai (ai + (d-2)/ii) 



(A.4) 



h'o 



2(d-2)(d-3) 
-2 (d - 2) (ao - Co) + 



3k^ai{ai-hi) (d - 1)^ a? + (d - 2) /i^ 



4(d-3) 



2(d 



, k'^ai{ai -hi) (d - 1)^ a? 

Co = 2 (ao - Co) - TTTi — 7:m — ^ + 



2(d-2)(d-3) (d-2) 



2 ' 



(A.5) 



where both the functions and their derivatives are evaluated at r = ro. In addition, there 
is a constraint 

2 (d - 2) / ao' + (d - 2) Co' [2 (d - 3) / + r /'] - 
2(d-3) (d-2) (ao-co)/' 



2 {d-2) f + rf 



+ {2 {d-2) f + rf) ho' + 



/2 



+ 



{d-l)k^r {d-3){d-l)\f 

2 {d-2) ~ {d-2) {2 {d-2) f + rf) 



Ol 



2 (d-l)r/a ; 
2 (d-2) 



/2 



0. (A.6) 



Second harmonic 
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The relevant equations are 



„ (d-2)/ + rf 

/ a2 H a2 + 

r 



-4A;^ + 



2 , 2(d-l)(d-3)/ 



/2 



[2(d-2)/ + r/ 



/1 2 



02 = S'rCaa , 



r.rf— 2 



(A.7) 



The sources read 



Src. 



0.2 



id-l){d-3)ffa [ , 1, ,2 
[2(d-2)/ + r/': 
(d-l)2(d-3)r/'3 



2 + 2 ^1 ^ 



+ 



oi - 
3k^2f-rf') 



Src 



■h2 



{d-2)[2 (d-2) f + rff 2 [2 (d-2) / + r/'] 
2(d-l)(d-3)//'(2/-r/')a;ai _^ ^,,2 



«i, 



3 (d - 2) r /ii ai + 



[2(d-2)/ + r/f 
4kHd-l)rfa[ , 4(ci-3)2(d-l)2//'' 2 



;i^"i + 



+ [2(d-2)/ + r/']"^+ [2{d-2)f + rf 

'e {A{d - 2)/ [(^2 - ^d + 5) f + {d - 4)r f] + {4d - 7^/ 



2 f/2\ 



[2{d-2)f + rfr 



af + 



2/-r/' 



These equations are subject to regularity boundary conditions at the horizon 



(Ai 



02 



2^2 
-2U-l-^^,a2 + 



(d-1)^ 



3fc2 



d - 2 2 (d - 3) 
, _ 8 (d - 2) 02 + ai [{Ad- 7)ai - 3 (d - 2) /ii] 2 



2 Sk^aihi 
"1+ 2(d-3)' 



2 (d - 2) (d - 3) 

where both the functions and their derivatives are computed at r = vq. 
The constraint is 

a^2(2/-r/0 2(ci-3) {d - I) f a2 

2 2{d-2)f + rf' (2{d-2)f + rff 

{d-l)rfa[^ 



(A.9) 



(d - 1) (d - 3) / ai a; 
(2(d-2)/ + r/')' "2(d-2) (2(d-2)/ + r/') 
(d-3) (d-l)V/'2 (d-l)fc2r 



+ 



(d-2) (2(d-2)/ + r/'f 2(d-2) (2 (d - 2) / + r /') 



of = 0. (A.IO) 



Mixed modes 

The mixed modes back-reaction appears in the T^ calculation starting from p = 2. 
The perturbations are governed by the equations 

{d-2)f + rf , 2 , 2{d-2)id-3)f'^ , 

f «ii + «ii - 2 A: an + — — — — —^ (an - cn) + Srca,, = 0, 

r [2 (d - 2) / + r /'] 
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/4 + M^M±lZ:,;,_2.^.„ 



2{d-3)f 



/2 



[2{d-2)f + rf'y 



(an -cii) + S'rccii = 0, 



/i^i'i H ^1 + 5'rc^ii = 0. 



(A.ll) 



where 

_ i{d - Iff'^ (-2(d - 3)(d - 2)f -id- l)rf'f + rV^2) 



+ 



{d-2)i2id-2)f + rf'y 
P (4((i - 2)/ (8(d - 2)2/^ - (5d^ - 27d + 33) rf/ - 5{d - 2)r2/'^)) 



(d-2)(2(d-2)/ + r/0=' 
(4d-7)r3f3 2, 4(d-l)//' 



a? 



a? + 



, 2/(2/-r/0 ,2 ^ 
7 OiOl + :7T^ ^r^r- nOi + 



+ 



(d - 2) (2(d - 2)/ + rff ' " 2(d - 2)/ + r/' "^"^ ^ 2(d - 2)/ + r/'^^ 
2/ (4(d - 2) (2^2 _ 5d + 1) /2 + (5^2 _ 6d - 7) r/'/ - (d - 3)rV'') /' ,, 



(2((i-2)/ + r/')' 

rff ,,2 , /(3rY2 + 2/(4/(d-2)2 + (4d-9)r/0) , 
f7 '^1 "I 77:71 „N ^ , 77T2 fli'^i) 



/i]^ai 



2{d-2)f + rf> 



(2(d-2)/ + r/') 



orCcj^j^ 



4(d - l)2/'2 (2(d - 3)(d - 2)/2 + (d - l)r/7 - r2/'2^ ^ 



(d-2)2(2(d-2)/ + r/0 



ai + 



+ 



2fe2 (2(d - 2)r3/'3 + (2 - d) (4(d - 2)p {{d^ + d - 7) f + {8 - d)rf'))) 



(8d - 21)r2// 



2ff>2 



{d-2f{2{d-2)f + rf'f 
2 4(d-l)//' 



-aid 



(d-2)2(2((i-2)/ + r/')' (d-2)(2(d-2)/ + r/ 

4/ ((d - 2) (5d2 - 14d + 5) /2 + 4(d - 2)dr/7 + r^ f'^) f 



2/2 



(d-2)(2(d-2)/ + r/')' 
,2 2/(r2/'2 + (d-2)/(2(2d-3)/ + 3r/0) , ^, 



h 



+ 



,'2 



Src 



toil 



2(d-2)/ + r/' 1 
2/ {rf - 2/) 
(d-2)(2(d-2)/ + r/')"^' 

2(d-2)cii(l+ ^•^ 



(d-2)(2(d-2)/ + r/'r 



ai/i'i + 



2{d-2)f + rf' 



k^ + 4aii 1 



(d-2)/ 



2{d-2)f + rf' 



e + 



2[(2d-3)r2/'2 + 4(d-2)/([d2-3d + 3] / + (d - 2) rf )] A;^ ^ 

H "T — _,,_,- :— r-Q Oi 



(d-2)(2(d-2)/ + r/'r 



+ 4 



(d-l)r/ 



m«i«i 



rf{rf'-2f) 



■ hi ai 



k\ 



Xd-2){2{d-2)f + Tf>)^'- [2(d-2)/ + r/f 
The horizon boundary conditions for these equations are given by 
, _ [(4d - 7)A:2 - 4(d - 3)(d - 1)^] af + 2(d - 2) [(fe^ - (d - 5)d - 6) an] 



(A.12) 



ni 



'11 



(d-3)(d-2) 
2 /2[(d-3)(d-l)2-(d-2)yt2]a2 



+ cn, 



d-3 



(d-2) 



+ (A:^ - d + 3) cii + an, 
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, _ 2fc^[(2ci-3)af + (d-2)(2an + (d-2)cn)] 
"^^1 " (d-3)(d-2) • ^^-^^^ 

where the functions and the derivatives are evaluated at r = rg. 
Finally, the constraint for the mixed modes reads 

, 2[{d-2)f + rf']a[, (d - 2) [2(d - 1)/ + rf] 4^ 2(d - 3)(rf - 2)f (an - cn) 
""^1+ 2{d-2)f + rf' ^ 2{d-2)f + rf' + (2(d - 2)/ + r/'f 

2(d-l)r/af 4((i-3)((i-l)/aia'i 2r//if 



(d-2)(2(d-2)/ + r/') ' (2{d-2)f + rff 2{d-2)f + rf' 

4 (d - 2)/ (-2(d - 3)/2 + (d - 5)r/7 + r^f'^) a^h{ 4((i - 3)((i - l)2r/'2a? 



{2{d - 2)/ + r/0' {d - 2) (2{d - 2)/ + rff 

B. Numerical issues and shooting parameters 



0(A.14) 



In this appendix we give some details regarding our numerical implementation and describe 
several tests that demonstrate its robustness. In addition, we summarize the numerical 
values of the shooting parameters in table 0. 

The ODEs encountered in this paper must be solved between the horizon and infinity, 
r G [ro,cxD). However, in practice we solve the equations only for r G [ro + ^iTmax]-, where 
e ^ ro and r^ax is some finite number. The first boundary is taken slightly outside tq be- 
cause the equations become singular at the horizon. To get the correct boundary conditions 
at this point we Taylor-expand our functions to the third order about the horizon and find 
the expansion coefficients analytically by plugging the expansion into the equations. This 
enables us to compute the functions and their derivatives at rmin = ^o + c- We checked that 
the meaningful numbers cited in this paper are independent of e (with better then 0.1% 
precision) provided e < 10~^. As for the outer boundary, its location is determined by the 
"empirical formula", Vmaxid) — 18 — (i in the range 4 < d < 14. Here too the meaningful 
numbers undergo less then 0.1% variations when we tried larger r^ax- 

In our shooting routines we integrate the equations from rmin , where one of the bound- 
ary conditions (b.c.) is used, toward rmax where we attempt to satisfy the asymptotic b.c. 
For k ^ this integration will in general diverge when r^nax -^ co since both the decaying 
and the growing solutions are present^''. Only for a specific value of the shooting param- 
eter the growing solution is eliminated, the integration converges and asymptotic b.c. is 
satisfied. Our method to filter out this growing solution follows that of Gubser Q: we 
use the integral of |x(r)p near the outer boundary, where x{r) is the mode we are solving 
for, as an indicator for how well the asymptotic b.c. is satisfied (for the true solution this 
integral must be negligible). The shooting and matching procedure is thus reduces to a 
minimization of that integral. 

The situation is different for the zero modes, as for them a finite solution exists re- 
gardless of the value of the shooting parameters. In this case we demand that the solution 



^For a finite rmax tiie result of this integration is a finite number tliat grows rouglily exponentially with 
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d 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


h,{l) 


.434 


.496 


.532 


.560 


.576 


.595 


.604 


.613 


.620 


.627 


.631 


ho{l) 


.796 


1.01 


1.24 


1.48 


1.73 


1.96 


2.21 


2.46 


2.71 


2.95 


3.21 


ao(l) 


-.232 


.015 


.506 


1.24 


2.23 


3.46 


4.95 


6.68 


8.68 


10.9 


13.4 


co(l) 


-.730 


-.921 


-1.14 


-1.38 


-1.62 


-1.86 


-2.11 


-2.36 


-2.60 


-2.85 


-3.10 


«2(1) 


-.40 


.690 


1.11 


1.66 


2.37 


3.22 


4.22 


5.43 


6.76 


8.26 


9.93 


h2{l) 


-.812 


-1.22 


-1.65 


-2.09 


-2.54 


-2.99 


-3.46 


-3.89 


-4.38 


-4.84 


-5.47 


aii(l) 


5.28 


7.74 


10.2 


12.5 


14.8 


17.1 


19.8 


22.2 


24.7 


27.1 


29.6 


cii(l) 


-2.96 


-2.92 


-2.90 


-2.88 


-2.86 


-2.84 


-2.88 


-2.89 


-2.89 


-2.90 


-2.90 


wn{l) 


2.96 


4.01 


5.04 


6.04 


7.03 


8.01 


9.08 


10.1 


11.0 


12.1 


13.1 



Table 7: The shooting parameters in various dimensions. We assume here that the normahzation 
for any p is same as for T^. The last 3 hues appear only for torii T^ with p> 2. 



vanishes asymptotically. This is easily accomplished for /iq, for whom "the second shot 
hits". However, as described in section |3.2| , the procedm'e should be iterated to make uq 
and Co vanish asymptotically. Moreover, because of the relatively slow power-law decay the 
location of the outer boundary at rmax^ given by the above empirical formula, is insufficient 
to ensure this asymptotic vanishing accurately enough. So, in practice, after integrating 
the equations ( [A.2| ) to r^ax we set the sources to zero (since they decay exponentially 
and are already very small) and continue to integrate to some large Vasymp (taken to be 



an order of magnitude larger than r„ 



In the vicinity of r, 



asymp 



we fit the solutions 



by expressions of the form Ooo + a/r'^~'^ and Coo + 7/r'^~'^, which is the correct analytic 
behavior in the asymptotic region [^]. Then we compute A = |coo| + laool, which is an 
effective measure for how well the boundary conditions are satisfied. To minimize A we 
iterate until A decreases below a certain tolerance. 

The numerical values of all shooting parameters appearing in this paper are listed in 
table 0. 

As soon as the solutions are obtained one must ensure that the constraints (p.^), 
( |A.6| ), ( [A. 10 ) and ( [A. 14 ) are satisfied by these solutions. In our case we verified that the 
constraints are satisfied with better then 0.1% precision. This together with other tests 
described above gives an idea of the numerical accuracy of our method. 
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